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Summary 


Three burns are covered: 760-1, 760-2, and 460-14. Temperatures, thermal 
radiation, etc. were analysed, but not air movements or pressure measurements. 
A standard empirical intensity-time curve was fitted to all of them, allowing 
intensity and time scale factors to be determined. The standard curve can be 
used in calculating the protection afforded by shelter and to calibrate the 
wood block system of instrumentation. It is hoped that the time scale factors 
constitute a measure of the burning rate and might ultimately serve to check 
various theories. However, discrepancies have been found between different 
elements and there were not sufficient data to resolve then. 


Introduction 


1. Flambeau is a project to study mass fires by burning arrays of piles of 
forestry fuel. All piles are ignited simultaneously by jellied petrol 
cartridges. 


2. This is the report of a hasty analysis of data received in the U.K. 

from the Riverside Laboratory, U.S. Forest Service, (1, 2 and 3) together with 
Countryman'’s article in Fire Technology (4). Three burns are covered: 760-1, 
760-2, and 460-14. Their vital statistics are shown in Table 1. 


3. The main point of the analysis was to get measures of the burning rate 

by fitting standard curves to the time variation of various elements:- Street 
temperatures, flame temperatures, thermal radiation, toxic gas concentrations, 
total and radiant fluxes, and weight loss rates. With further theoretical work 
it might prove possible to analyse soil temperatures in a similar way, but this 
has not yet been done. No use has been made of meteorological data, air move- 
ments or pressure measurements. The space variation of burning rate shouwid 
provide a check on various theories of interaction (5) and the standard curves 
should help in assessing the protection afforded by shelters. 


Method of Analysis 


4, In most cases the raw data were extremely irregular and variable. A 
preliminary smoothing was carried out; in the case of graphical data by 
pencilling in an average curve by eye, and in the case of numerical data by 
taking means over a period of a few minutes. 


5. The analysis was based on the assumption of theories 20a and 20b (5), that 
all these intensity-time histories can be represented by one or more standard 
curves. To do this, it is neceseary to fit:- 


(a) The time scale factor 

(b) The intensity scale factor 
and (c) The time scale origin 
(There is no difficulty about the intensity scale origin, because this is 
either zero, as with Weight loss Rate, or the ambient pre-burr value, as with 
temperature). The curves all seem to be skew, with a long tail at late tines, 
so a Pearson Type III distribution was tried. This did not fit satisfactorily, 
and it was decided to fit a wholiy empirical curve. 
6. The two scale factors can obviously be easily fitted by plotting on 
log:log paper, but this leaves the problem of the time scale origin. The obvious 
eolution is to take this ae the ignition time, but this makes the curve depend 
on the early stage of fire growth after ignition, which is not of much interest 
in Flambeau, and can be pretty variable. The next possibility is to take it as 
the time of peak intensity, but this has 4 disadvantages:- 


(a) The curve is split into two separate segments, with the peak ituelf 
off the paper. 


(b) Many of the curves finish before the peak is reached 
(c) Some havevery flat peaks 


and (d) Some have double peaks. 


mae 


The solution adopted was to measure time from the moment when the intensity 
first reached 2 of the peak intensity (eee Figure 1). This is not an ideal © 
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solution, but it ie believed that it should reduce such errors to acceptably 
low levels. For example, so long as the curve goes on for a substantial time, 
even if it does not reach the peak, it should be possible to guess the peak 
inteneity with sufficient accuracy that there cannot be any great error in the 
time of 2 peak intensity. 


7. As stated above, Ip was plotted against t ~ t(ZIp) on log:1log paper. The sls 
factors were characterised by Ip, the peak intensity, and Ot (3 I p), the 
duration at more than 3 the peak intensity - written as A. 


8. Out of a total of 49 curves two have been rejected. The first was the 
flame temperature for Pile B3, Burn 760-1, rejected because the >eadings 
continued for such a short time. The second was the flame temperature for Pile 
C10 (or Position 29), Burn 460-14. This was based on one seneor only, the 
results were quite inconsistent with flux measurements on the same piie, and 
most of the time when the pile was burning fiercely the sensor was merely 
measuring the temperature of cold air (2). 


9. There is no great difficulty about fitting if the curve has a clear and 
satisfactory peak. Im the case of curves that have not, a further aasunption 
is needed to be able to fit; but afterwards these aleo can contribute to the 
shape of the curwe. Because of this feedback, fitting has beea done by a 
process of successive approximation, as shown in Figure 2. 
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10. There are a number of different further assumptions that can be used:- 


(a) That 4 ts the same for al) curves at the same point, or relating to 
the sa®e pile. This is justified by the observation that in such 
cases the peaks are normally at the same tine. 


(>) That the time integral of intensity is the same at all piles, whether 
they burn fast or slow. This is more or less equivalent to the 
assumption that the total energy release is the eame for all piles, 
and hence that the efficiency of combustion is the eame. The integral 
is equal to 1.68 4.Ip for the curve which has been fitted, so the 
value of A .Ip must be the came for the same type of observation. 
This assumption has been used, in particular, to analyse street 
temperatures for 460-14 and flame temperatures for 760-1. 


(c) That the ratio of street temperature integrals at two heights is the 
same everywhere within a given burn. 


(d) That the ratio of the integralsof two different effects is the same on 
different burns. This assumption seems fairly doubtful (see Paragraph 
17 below), but may provide a certain degree of guidance. 


(e) The ratio between street temperature integrals in different arrays can 
be calculated by the inverse square law - theories 3a and 3c (5). 


This give the ratios 760-2 and 460-14 
om = 7.19 Feo. * 10.01. 


When several of these methods were available a weighted mean was taken, the 
weight depending on the number of observations and the directness of the 
inference. 


Results 


11. Figure 3 shows the curve which has been fitted. The dots indicate the 
degree of spread around it, which is also shown by the lines of + 1 and 2 
Standard Deviations. The curve has also been tabulated in Table 2. A is equal 
to 2.80 x peak tise (measuring from time of ~ peak), and, as stated, the 
integral from 0 to invinity is 1.68 4 . Ip. However, there is a shortage of 
points towards the end of the curve, and depending on the slope chosen, the 
integral could be almost indefinitely larger (it could not be very much cerlier). 
A single curve has been used to fit all the observations, i.e. theory 20a (5) 
has been followed - but see the end of Psragraph 17. 


12. Table 3 shows the average height variation of the integrais of street 
temperatures. 


Table 3 
Height Variation of Street Temperatures 


460-14 - integral | 19,000] 25, 300] 30, 100/32, 800/34, 400! 35,400 


| 

" - normalised to 100@7'| 75/ 100/ 119] 130) 136} 140 
760-2 - integral } 2 > |21.500 | 
" o" - normalised to 100 @ a - - | 67 | 


In this table, the figures for 460-14 are based on 20 different curves and sre 
probably moderately reliable, while those for 760-2 are based on only ome curve 
at each height and are therefore quite unconfirmed. 


13, Table 4 shows, for each of the 47 remaining curves, the values of 4, Ip 
and the time integral, calculated by the methods of Paragraph 10, In addition, 
these values are shown calculated by direct fitting, as in Paragraph 9, where 
this was possible. These, fitted only from the curve itself, are more relevant 
when checking consistency and the assumptions of Paragraph 10. The times of the 
more gee peaks (measured from ignition time) and the gap between ignition 


and t (4 Ip) are aleo quoted. 


14, The duration O should be inversely proportional to the burning rate, 60 
it seems reasonable to use this as a measure of it. Table 5 shows the mean 
duration at each location, where this is known. 


Discussion 


15. Figure 4 shows the durations for 760-1. As would be expected with this 
wide separation, there seems no evidence of any important variation or inter- 
action. The dry-bulb temperature results are of considerable interest. The 
upwind station, P3, has somewhat lower integrale than P1 and P4, which are at 
the same distance from the array but downwind or crosswind. Moreover, the 
interior station P2 has an integral 2.68 times the mean of P1, 3 and 4 (using 
the directly fitted figures), which compares with a ratio of 2.96 calculated 
on the assumptions of Theory 3 (5). This agreement tends to confirm these 
assumptions, viz:- 


(a) The temperature rise at a point is the sum of the contributions from 
each of the piles. 


(>) The contribution is proportional to the burning rate (assumed equal 
for all 9 piles). 


(c) The contribution is inversely proportional to the square of the 
distance (confirmed for 80' ~ 450"). 


16. Figure 5 shows the Mean Durations for 460-14. A set of contours have been 
added based entirely on the street temperatures; if the other measurements are 
for the moment ignored it can be seen that:- 


(a) There are important variations from place to place. 


(b) On these 6 points, the spacial variation seems very sooth and 
regular, rather than random variation between piles; but the other 
measurements tend to contradict this. Om the other hand, qualitative 
accounts by observers support it. 


(c) The variation is not radial and does not seem to be due to wind 
effects - the ambient wind was alwost a dead cals. 


Hence interaction at peak whether due to wind (Theory 4 (5)) to gas composition 
and temperature (Theory 5) or to thermal radiation (Theory 6) does not seem to 
explain the results, and yet there did seem to be interaction, from (b). The 
explanation might possibly be that interaction was by ignition or re-ignitior, 
before peak (Theory 2c). 


17. The 3 remaining Durations for 460-14 are all relstively low, but do not 
agree with the street temperatures or with each other. Possible explanations 
include:- 


(a) Despite the evidence in Paragraph 16, there are large random 
variations from pile to pile. 


(>) The weight loss measurements, which were only made on z of one pile 
and do not give a very good shaped curve, are to be discounted. 
However, the integral 6700 lbs. was very close to the total weight 
on the platfora, 7000 lbs. 


(c) Similarly the flux measurements might be discounted. 


(d) The Thermal Radiation Duration was low because the edge of the array 
burnt faster. This would suggest some version of Theory &, or 
Theory 5 if the effect of fire gases is to depress the burning rate. 


(e) The Thermal Radiation Duration was low because smoke accumulated during 
the burn, greatly reducing the intensity at a later stage and moving 
the peak to an earlier time. 


(2) The Thermal Radiation Duration was low because early on there was a 
big flame and intense sideways radiation; but later, although the 
burning was still continuing at a great rate most of the heat was 
emitted by upward radiation and by convection - the radiometer could 
see little but the cold outsides of logs. 


(g) The street temperature durations were higa because hot gas was steadily 
accusulating in the streets during the burn. 


(bh) The street temperature durations were high because until the low 
level air had been used up hot air did not come down into the streets 
(Theory 5a). Against thie, there does not seem to be any evidence 
that the temperature rose faster at 50’ than at 3.5'. 


It does not seem possibile to distinguish between these explanations on the 
evidence now available. It should be noted that the last 4 involve systematic 
differences between the elements - i.e. Theory 2Ob rather than 20a. 


18. There are not many durations available for burn 760-2. It seems that this 
may have burnt eomwhat faster than 760-1, and much faster than 460-14. As 
before, the duration for radiation is a lot shorter than for street tempera- 
tures, which tends to confirm that there is a real effect; but it might be either 
a radiation/temperature effect, or an outside/inside effect (the radiomcter was 
on the outside of the array, whereas street temperatures were taken inside, in 
both cases). In this burn the variation of toxic gases (over a pile) was not 
inconsistent with that of flame temperatures, whereas in 460-14 the toxic gases 
(in the street) did not sees to correlate well with street temperatures. As 
shown in Table 3, the street temperatures seemed to fall with height instead 
of rising as in 460-14, If this observation is valid, it suggests that 
vertical circulation was much more important in this case ~ leading perhaps to 
some version of Theory 5. 


19. Table 4 also shows data for 5 curves from Countryman's paper (4). These 
have not been used in the analysis, since it is not known where they were, or 
in which burn, but they were all good shaped curves. It is therefore interest- 
ing that the standard curve, which was derived without reference to them, 7its 
them all extremely well. 


Recommendations for Future Work 


20. Not enough observations were available here to answer many of the vital 
questions. In future, it is suggested:- 


(a) We need an extensive set of measurements all on the sase element 
throughout the array. This could be provided by the wood biock 
system, which could be calibrated using the stancard curve. 


(b) We also need timed readings at a limited number of points. It is 
very important that these should go on for a long time after the 
peak. 


(c) The discrepancy between temperature, thermal radiation etc. needs 
clearing up. This could be done by intensive observations on singe 
isolated piles. Thermal radiation usually gives very good curves - 
if only we can rely on it. Weight loss is very fundamental - if it 
can be got to provide good curves. 


(da) To check interaction effects, one or wore isolated piles should 
always be burnt at the eame time as any array. 


(e) We need to know whether a small proportion of fire gases in the input 
air increases or decreases the burning rate. This couid be discovered 
in the laboratory. 
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Table 1 


: Particulars of the three burns 


Pile dimensions 


Pile weight 


Pile energy (BTU) 


Width of streets 
(pile separation) 


Number of piles 


Size of array 


Ignition 
Date/tine 


lapse rate 


Air temperature 


Relative humidity 


wind speed 


Wind direction 


Ground 


Moisture content 
~ fine fuel 


Moisture content 
- heavy fuel 


Rurning rate 


Burning efficiency 


47'x &7'x 7! 


20 tons 


340 x 10° 


415" 
pee Shee) 
DA XeoTAs 


Jan. 31,1964/1005 


39°F 
2m 
3-5 knots 


Easterly 


Dry 


LPtx Lite Ft 
20 tons 


340 x 10° 


25u 


AGUS 1B = 2k 


1271 "x 1271" 


Dec.6,1965/1215 


Stable-inversion 
-unstable 


13°F 
61% 
0 - 3 knots 
Variable 


Under 4"-6"snow 


1b 


47S TS a 


20 tons 


340 x 10° 


25) 
6x6 = 36 
407'x 407! 


May 15,1964/0830 


Unstable 


Table 2; The stander intensity-time curve 


I (arbitrary units) | 


0.71 


t (mins) |I (arbitrary units) 


(7 


0.8 


0.65 
0.59 
0.52 
0.44 
0.38 
0.33 
0.29 
0.26 
0.21 
0.166 
0.130 
0.103 
0.080 


0.067 


0.045 
(0.017) 


(0.004) 


Notes 1. These values have been calculated for G= Ip=1. For other 
values of 4 and Ip, scale proportionately. 


2. All times are measured from t(3Ip), the time of + Peak Intensity. 


3. Values from t = 6 onwards are extrapolated and highly uncertein. 
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(2) 411 temperatures in P, Pluxes in BTU/ft2/seo, concentrations in %, Thermal radiation in BIU/it/nr. 


78+ means the ou-ve 


Table 5 _: Mean durationsat various locations 


Duration A (winutes} 


Flame temperature 


Stetions Pi,2,3,4 Dry-bulb temperature 
Pile ci0 Total & radiant flux 

(DAO Weight loss rate 

R - 50' from array Therwgal radiation 


Station T6 Street temperature 


Wi ~ 50‘ from array Thermal rediation 


Centre Street temperature 


-10- 


-bl- 


Intensity 
(arbitrary units) 


Q.1 


Cf lS ee eee 


0-01 


Time (ninutes) 
Fig.3 THE STANDARD INTENSITY-TIWE CURVE, SHOWING THE SCATTER OF POINTS ABOUT IT. 
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Fig. 4 Burn 760-1 


DURATIONS AT VARIOUS PLACES 
8 


1 
ony ; | ee i ‘a 


iat 
[ 
[| 


P, @20 


o1 De 


Burn 460-14 


Fig. 5 


DURATIONS AT VARIOUS PLACES 


a) 
O 
fa} 
a) 
0 
O 
a) 
QO 
0 
o 


BA 8/8 00> 6o0, eo eeete re ese see 


i 


11 


iS) (a) ey (eh ja} je) 
2) (o) tet 8) oa) 1G. 
je! jl je) oe) 
vale ear al (oi 
le) (aire) se) el 
ee) (oath (cl) 
ala) (e) eo) ie) (a, 
[B) (a) 45), ae) 1) 
(a) so) fe) aie) ea) so 
ey o)e(a) eal (a) (a) 
Oa, eee 
BF OSes ers 
yee el 
ee) feel oO) aay (o) 
Bb Bea 2. eB 
en OJ J a ap 
O° Oe) bare 
ogoagaagaocada 
pate Ss ee 


— 


Od Baers Be Rea are eee 


a 
4 


18 


Note: 


Dots (@) are street temperature measurements, crosses (x) are others. 


Contours based only on street temperatures. 


